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SECTION  1.  OVERVIEW 


This  system  allows  the  collection  and  display  of  neurophysiological  data  during  the 
time  the  experiment  is  being  run.  The  system  was  initially  designed  for  collecting 
both  intracellular  and  extracellular  evoked-response  information.  A  25-ms  sample 
from  two  channels  of  analog  information  is  taken  periodically  and  processed.  The  two 
channels  of  information  are  sampled  as  closely  to  simultaneously  as  the  instruction 
timing  will  allow  (using  an  LSI-11,  30  ys).  The  operating  mode  and  intersample  period 
may  be  controlled  by  the  computer  through  a  user-supplied  time  interval.  It  also  may 
be  driven  by  a  Schmitt  trigger  through  timing  pulses  during  an  experiment  or  from  a 
recording. 

The  program  runs  in  two  parts.  The  first  part  (SETUP  portion)  allows  the  user  to 
define  areas  of  interest  on  the  sampled  waveform.  The  second  part  (experiment  RUN 
portion)  collects  samples,  looks  at  the  areas  defined  in  the  setup,  does  some  prelimi¬ 
nary  analysis,  and  displays  the  data  on-the-fly  at  the  terminal.  The  setup  portion 
collects  samples  and  displays  the  average  waveform  for  both  channels  simultaneously 
on  the  CRT.  The  user  may  then  specify  time  windows  (which  bracket  features  of 
interest  on  the  waveform)  by  setting  cursors  on  the  waveforms  displayed. 

The  user  has  control  over  the  following  parameters: 

1.  the  operating  mode,  either  event-driven  or  computer-driven 

2.  the  number  of  signals  averaged  during  the  setup  portion;  the  time 
domain  amplification  shown  on  the  screen  for  setup;  defining  areas  of 
interest  on  the  waveforms,  using  time  windows 

3.  determining  whether  voltage  measurements  will  be  made  relative  to 
zero  levels  or  as  absolute  differences  within  the  defined  time  windows 
(zero  levels  are  continuously  updated  to  negate  the  effects  of  DC 
drift) 

4.  scaling  on  the  output 

5.  saving  the  results  by  supplying  a  file  name  for  the  data  to  be  stored 
when  asked 

All  programs  are  designed  in  a  strictly  modular  fashion,  which  allows  the  quick  and 
easy  addition  of  analysis  modules  as  required  by  the  user.  Modules  are  written  in 
either  FORTRAN  IV  or  MACRO  (PDP-11  assembler).  Certain  system  routines  are 
used  from  the  RT-11  operating  system;  thus  the  use  of  another  operating  system 
requires  the  replacement  of  these  calls  with  similar  ones  from  that  system. 

The  programs  use  named  common  areas  in  memory  to  allow  new  modules  to  access 
only  those  variables  and  flags  that  are  needed.  This  technique  effectively  "hides" 
other  variables  and  minimizes  side  effects. 

Precision  in  the  data  handling  is  as  follows.  An  input  signal  ranging  from  -5  V  to  +5  V 
is  passed  into  a  12-bit  A/D  converter.  This  results  in  a  discrimination  of  0.024% 
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provided  that  the  full  10-V  range  is  used.  The  data  stored  in  the  data  files  reflect  this 
precision.  The  average  waveform  plotting  is  correct  to  1%  of  the  maximum  range  on 
input  due  to  scaling.  The  data  graphing  during  the  experimental  run  is  correct  to  2% 
of  the  range  specified  by  the  user  due  to  scaling. 


HOW  TO  USE  THIS  DOCUMENT 
General  User 


A  potential  user  should  review  the  Figures  and  Output  section  along  with  this  overview 
to  determine  applicability  of  this  package  to  his/her  needs.  The  user  should  also 
consult  the  hardware  and  software-module  description  (see  Equipment  section  below). 
When  the  package  is  assembled  and  ready  to  run,  the  procedure  description  portion 
(Section  2)  should  be  kept  nearby  for  reference. 

Programmer 


All  the  code  necessary  to  run  this  task  is  included  in  this  document  in  the  appendices. 
The  internal  documentation  of  the  code  and  the  modularity  of  the  design  should  facili¬ 
tate  adding  to  the  main  code.  If  more  off-line  analysis  routines  are  needed,  consult 
Appendix  E  for  the  data  file  structure. 


EQUIPMENT 

Equipment  needed  to  use  the  system  is  listed  in  the  following. 
Hardware 


Hardware  equipment  needed  is  CPU,  video  terminal  with  graphics  capability,  A/D 
device,  Schmitt  trigger  device,  programmable  clock,  signal  source  ±5  V,  and  parallel 
output  port. 

Software 


Software  equipment  needed  is  the  RT11  operating  system  and  other  program  modules 
(RT11  is  a  trademark  of  the  Digital  Equipment  Corp.): 

swante.for 

swanse.for 

swanru.for 

colect.mac 

getbuf.mac  (supplied  by  AFRRI) 
stimul.mac 

plot55.mac  (supplied  by  DEC) 
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CHECKLIST  FOR  PROPER  SETUP  OF  SYSTEM 

1.  Ensure  all  equipment  is  connected  as  desired. 

2.  Turn  on  the  LSI,  terminal,  disc  drive,  etc. 

3.  Boot  system, 

4.  Turn  on  line  clock. 

5.  Build  task  if  not  done  yet. 

6.  Start  task  (RUN  DKlrALLFOR).  The  program  is  up  and  running  if  the  CRT  dis¬ 
plays  Figure  1. 


WILL  THIS  EX.  BE  SCHMITTTRIGGER  DRIVEN?  (Y  OR  N) 
Figure  1.  First  question  asked  when  program  starts 

7.  Answer  all  prompts. 
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SECTION  2.  PROCEDURES  FOR  OPERATION 


SETUP  PROCEDURES 
Baseline  Procedure 


The  following  steps  must  be  performed  to  calibrate  the  setup.  If  the  amplifier  gain  is 

changed  after  these  steps  are  performed,  the  results  will  not  be  accurate  in  magnitude 

(although  still  correct  in  direction). 

1.  When  the  line  is  at  0  mV  on  the  first  channel,  hit  <cr>  (carriage  return). 

2.  When  the  line  is  at  5  mV  on  the  first  channel,  hit  <cr>. 


NOTE:  The  <cr>  must  be  hit  while  the  line  is  at  5  mV;  therefore, 

the  longer  the  standard  pulse,  the  better. 


3.  Repeat  steps  1  and  2  for  the  second  channel. 

Amplification  Procedure 

After  drawing  the  axes,  the  program  will  wait  for  you,  the  user,  to  tell  it  how  much 
time  to  plot  in  the  screen  width.  An  amplification  of  1  will  plot  25.6  ms,  an  amplifi¬ 
cation  of  2  will  plot  12.8  ms,  an  amplification  of  3  will  plot  8.5  ms,  etc.  (i.e.,  25.6 
ms/amplification  factor). 


You  may  change  the  amplification  as  much  as  you  desire  by  entering  another  value 
each  time  the  amplification  prompt  comes  up.  If  you  hit  a  <cr>  without  a  value,  the 
program  will  then  proceed  to  the  cursor-placing  routine. 

Cursor-Placing  Procedure 

Cursors  must  be  placed  on  the  graphs  to  delimit  time  windows.  One  pair  comprises  a 
time  window  within  which  the  analysis  will  take  place. 

Cursors  are  placed  in  the  following  way: 

1.  Choose  a  graph  as  directed. 

2.  Enter  a  number  (1-256). 

3.  Hit  <cr>.  A  cursor  will  now  appear  at  the  position  requested. 

4.  Do  you  want  to  save  this  position  as  a  time  window  delimiter? 

If  yes,  hit  <cr>  and  then  go  to  step  2. 

If  no,  go  directly  to  step  2. 
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5.  When  you  are  through  placing  cursors  and  defining  time  windows,  enter  a  2  and  hit 
<cr>.  The  cursor-placing  procedure  is  now  complete. 


RUN  PROCEDURES 
Scaling  Queries 

The  maximum  values  must  be  supplied  as  asked,  including  the  units.  This  procedure 
determines  the  scaling  factors  for  the  output.  The  number  of  trials  desired  must  be 
supplied  as  asked.  The  program  will  now  run  on  its  own,  either  driving  the  prep  with 
the  stimulator  or  through  a  Schmitt  trigger  arrangement.  The  program  will  auto¬ 
matically  make  copies  of  the  data  on  the  screen  as  necessary. 

After  all  trials  have  been  run,  the  program  will  prompt  for  a  name  for  a  file  to  store 
the  data.  Once  this  is  done  and  the  data  are  stored,  the  last  question  the  program  asks 
is  "Do  you  wish  to  run  more  trials  with  the  same  set-up?"  If  the  answer  is  yes,  the 
program  will  loop  back  into  the  "run"  mode  again  without  going  through  a  setup  pro¬ 
cedure.  This  feature  can  be  useful  in  determining  those  scaling  values  that  give  a  good 
display  of  the  data.  If  the  answer  is  no,  the  program  will  end  and  display  all  the  data 
files  that  are  stored.  At  this  point,  other  analysis  routines  can  be  used  on  the  stored 
data.  (Please  note  that  Figures  11  and  12  show  a  sample  output.) 
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SECTION  3.  FIGURES  AND  OUTPUT 


This  section  contains  examples  of  the  questions  asked  to  establish  parameters  and  the 
valid  responses  to  those  questions.  Also  the  type  of  output  and  displays  are  shown. 
The  figures  included  show  actual  output  as  seen  on  the  CRT. 


SETTING  INITIAL  PARAMETERS 

Figure  2  lists  all  initial  questions  that  establish  the  setup  parameters.  For  acceptable 
answers  to  the  questions,  refer  to  Table  1. 


WILL  THIS  EX.  BE  SCHMITT  TRIGGER  ORIVEN?  (Y  OR  N)  N 
ENTER  2  CHANNELS  FOR  INPUT  .1.2 
WHEN  CHANNEL  1  IS  ZERDED,  HIT  (CR) 

1  9B7 

WHEN  5MV  STANDARD  APPLIED,  HIT  (CR) 

1  929 

WHEN  CHANNEL  2  IS  ZERDED,  HIT  (CR) 

2  1045 

WHEN  5  MV  STANDARD  APPLIED,  HIT  (CR) 

2  877 

#  OF  SIGNALS  TD  BE  AVERAGED:  3 
DD  YOU  WISH  TO  DUMP  THE  BUFFER?  (Y  DR  N)  :  N 
ENTER  THE  TIME  BETWEEN  SAMPLES  (SEC)  :  5 


Figure  2.  Setting  initial  parameters 


Table  1.  Acceptable  Answers  to  Questions  in  Figure  2 


Question  or  Requirement 

Acceptable  Answer 

Number  of  signals  to  average 

Any  number-^  3 

Enter  two  channels  for  input. 

Any  two  channels  0-14  separated  by 
comma.  First  position  refers  to 
stimulating  presynaptic  channel  and 
upper  graph.  Second  position  is 
postsynaptic  channel.  Be  sure  chan¬ 
nels  selected  are  ones  signals  are 
going  into. 

Do  you  wish  to  dump  buffer? 

This  option  allows  user  to  test  A/D 
converter.  If  chosen,  last  buffer  of 
digitized  data  is  output  to  CRT. 
Program  execution  then  halts. 

Enter  time  between  samples. 

Enter  time  as  integer.  This  question 
will  not  appear  if  Schmitt  trigger 
driven. 
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Each  time  a  sample  is  taken  during  the  setup  portion,  the  word  "WORKING"  will 
appear  on  the  screen  (see  Figure  3).  The  number  of  times  it  appears  corresponds  to 
the  number  of  samples  to  average.  The  time  between  appearances  equals  the  time 
between  samples. 


#0F  SIGNALS  TO  BE  AVERAGEO  :5 

00  YOU  WISH  TO  OUMPTHE  BUFFER? 

(Y  OR  N)  : 
WORKING 

N 

1 

WORKING 

2 

WORKING 

3 

WORKING 

4 

Figure  3.  Output  seen  on  CRT  while  sampling 
during  setup  portion 


After  sampling,  the  graph  will  be  drawn  without  plotting  (Figure  4).  It  requires  an 
"amp  factor," 


NOTE:  EACH  POINT  =  0.1  MS 

CURSOR  X  j 

1  { 

2  j 


CURSOR  X 

1 

2 

3 

4 


AMP  FACTOR. 


CURRENT  POSITION 
Y  = 

SE.  Y  = 


Figure  4.  Machine  prompts  user  for  an  "amp  factor." 
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USING  AMP  FACTOR  OPTION 


The  "amp  factor"  option  allows  the  user  to  expand  and  amplify  the  time  domain 
plotted  on  the  screen.  Figures  5  through  8  demonstrate  the  effect  of  various  choices. 


Figure  5.  Sinusoidal  waveform  at  amp  factor  =  1.  Compare  this  plot  with 
Figures  6,  7,  and  8  to  see  how  this  option  works. 


Figure  6.  Sinusoidal  waveform  plot  at  amp  factor  =  2.  Compare  this  with 
Figures  5,  7,  and  8. 
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Figure  7.  Sinusoidal  waveform  plot  at  amp  factor  =  5.  Compare  this  with 
Figures  5,  6,  and  8. 


Figure  8.  Sinusoidal  waveform  plot  at  amp  factor  =  10.  Compare  this  with 
Figures  5,  6,  and  7. 
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Now  the  graph  on  which  to  place  cursors  can  be  chosen  (see  Figure  9). 


Figure  9.  User  ready  to  define  time  intervals 


The  cursors  have  been  placed  on  both  graphs,  and  data  are  ready  to  be  collected  (see 
Figure  10). 


Figure  10.  User  defining  time  windows 
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As  the  experiment  is  run,  the  results  are  plotted  on  the  CRT  as  shown  in  Figures  11 
and  12. 


P  LATENCY  |  TOTAL  ELAPSEO  TIME: 

0  STIMULUS— )MAX.  1 - .1  O' 7:0 

S  NEG.EXCURS  1 1- :  j 

T  (MSEC)  I  j 

s . 1 . ; . 

Y  OELTA  VOLTAGE  t  +  +  +  +  +  +  +  +  +100+  +  +  +  +  +  +  +TRIALS 

N  MAX.  :  : 

P  NEG.EXCURS  10- :  j 

S  (MV)  !  j 


. + . + . + . . . . + . +“*V(5*0  + . + . . ’+fmALS 

P  OELTA  VOLTAGE 
R  CRSRS  1-)2  125- 

E  (UV) 


S  OELTA  VOLTAGE 
Y  CRSRS  3-)4 
N  (MV) 

P 

s . 


'+  +  +  +  +  100  +  +  +  +  +  +  +  +  TRIALS 


+  +  +  +  +  100  +  +  +  +  +  +  +  +TRIALS 


TRIAL:  30 


Figure  11.  Sample  output 


Figure  12. 


Sample  output 
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APPENDIX  A.  BUILDING  AND  LINKING 


BUILDING  THE  TASK 

The  following  sequence  is  used  to  compile  and  link  the  task: 

Compile 

FORTRAN 

FOR/EXTEND/OBJ:DKl:ALLFOR/LIST:DKl:  DK1:(SWANTE+SWANSE+SWANRU) 

This  causes  all  the  FORTRAN  programs  to  be  compiled  into  an  object  module  called 
ALLFOR.  This  module  is  placed  on  DK1:;  it  is  assumed  that  the  FORTRAN  source 
code  is  also  on  DK1:. 

macro: 

MAC/OBJ:DKl:  DK1:GETBUF/LIST:DK1: 

M AC/OBJ: DK1:CL  DK1:C0LECT/LIST:DK1: 

MAC/OBJ:DKl:S  DK1:STIMUL/LIST:DK1: 

M AC/OBJ: DK1:PL  DK1:PL0T55 

The  macro  programs  are  compiled  separately,  and  in  some  cases  the  names  of  the 
object  code  files  are  shortened  for  convenience. 


LINKING 

LINK/EXE:DK1:/MAP;DK1:  DKl:(ALLFOR,CL,PL,GETBUF,S) 

It  is  usually  convenient  to  create  an  indirect  command  file  that  performs  the  tasks 
outlined  above. 
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APPENDIX  B.  EXAMPLE  OF  HARDWARE 


The  hardware  used  for  the  development  of  this  system  is: 


CPU 

-  LSI-11 

Terminal 

-  VT55-FA 

Schmitt  trigger/clock 

-  KWV-11A 

Mass  storage 

-  RK05 

A/D 

-  ADV11-A 

Parallel  output  port 

-  DRV  11 
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APPENDIX  C.  LISTINGS  FOR  FORTRAN  PROGRAM 


C.1  SWANTE .FOR 


REAL-TIME  INTRA/INTER  CELLULAR  EVOKED  RESPONSE  DATA  COLLECTION 
PROGRAM.  MAY  BE  SCHMITT  TRIGGER  DRIVEN  OR  IT  MAY  DRIVE  THE 
PREP  ACCORDING  TO  A  TIME  PARAMETER  WHICH  IS  PROMPTED  FOR. 


MASTER  PROGRAM 
0  CALL  SETUP 

CALL  PLOT55C  9,0,0) 

CALL  PLOT55( 10, , ) 

CALL  PLOT55 (9,35,11) 

CALL  PL0T55O2,  ,  ’END  OF  SET  UP;HIT  <CR>  TO  CONTINUE’) 
CALL  PLOT55 (9,35,14) 

CALL  PLOT55( 12, , ’ENTER  S  TO  RESTART’) 

CALL  PLOT55(2,512, ) 

ACCEPT  1000,  S 
IF  (S.EQ. 'S’ )  GO  TO  10 
CALL  RUNIT 
000  FORMAT(AI) 

END 


C . 2  SWANSE .FOR  -  SET-UP  PORTION 
c 

c  25  march  80 

c 

c  routine:  setup 

c 

c  routine  for  the  initial  set-up  of  the  experiments 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


this  routine  receives  an  initial  5mv  square  pulse 
and  calculates  the  function  that  maps  the  real 
voltage  across  the  membranes  from  the  digitized 
input,  it  then  uses  this  value  later  to  output 
the  potential. 

the  parameters  passed  include  the  slope  ’m'  and 
the  y  intercept  ’b’  of  this  function. 

following  this  action,  the  program  receives  and 
averages  'n  ’  signals  and  presents  these  vs  time, 
the  experimenter  then  places  four  cursors  along 
the  signal  trace  to  indicate  the  points  between 
which  he  wishes  the  dv  measured. 

main  prog 

/ 

setup  runit 

/  ~  ! 

stdsig  sigavg  csrplc 

/  /  / 

colect  getbuf  axes  trcplt  store  disp 

/ 

scale  plot 
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c 


subroutine  setup 


Appendix  C 


C 

call  plot 55(2, 1+2 +4+32+6 4+ 128+ 25 6+512,)  lenable  graphics 

call  plot55( 13,72, )  !  clear 

call  plot55(  13,74,  )  !  screen 

call  stdsig(l)  !  find  mapping  function 

call  stdsi g(2 ) 

call  plot55(  13,72,  ) 

call  plot 55 (  13,74,  ) 

call  sigavg  liampl  is  defined  in  trcplt 

call  csrplc 

call  copy  !copy  screen  for  later  reference 

c 

c  set  up  is  now  complete 

retu  rn 
end 
c 

c _ 

c _ _ _ 

C  "  . 

c  routine  for  using  a  standard  signal  to  generate 

c  the  mapping  function  for  the  digitized  values  to 

c  their  mv  equiv.(for  output  later) 

c 

subroutine  stdsig(i  ) 

common  /scalar/ic(2,4) , ib(2),m(2) 

common  /schmi t/ist , izzz 

common  ichtb 

integer  volt0,volt2 

byte  ichtab(2) 

equivalence(ichtab, ichtb) 

if  (i.eq.2)  go  to  20  Iskip  if  already  answered 
10  wr i te  ( 5 , 1 000 ) 

1000  format ( 1 $will  this  ex.  be  schmitt  trigger  driven?  <y  or  n>  :',$) 
accept  1500, ist 
1  500  format  ( 1  a  1 ) 

i f ( i st . ne . 'y ' . and . ist . ne .  ' n  1  )go  to  10 
if  (i st . eq .  *y 1 )  ist= 1 
if  (ist . eq .  'n  ’  )  i st=0 
wri te  ( 5 , 2000 ) 

2000  format ( 1 $enter  2  channels  for  input  :') 

accept  2500,  ichtab ( 1 ) , ichtab ( 2 ) 

2500  format ( 2i 3 ) 

20  write  (5,3000)i 

3000  format  ('$when  channe 1  1 , i 3  ,  '  is  zeroed,  hit  <cr>'  ) 
accept  3500 
3500  format  (al ) 

j  =i chtab (i  ) 

call  colect ( voltO , j )  !get  sample  from  a/d  @  0  mv 
write(5,4000)j, voltO 
4000  format(2i10) 

wri te (5 , 4500) 

4500  format  ('$when  5mv  standard  applied,  hit  <cr>' ) 

accept  3500 

call  colect (volt2 , j )  ! get  sample  from  a/d  @  5  mv 

wri te(5 ,4000)j, volt2 
do  30  ii=1 ,32000 

30  continue  ! pause  for  a  moment 

itry=.2  *  ( volt2- voltO ) 
if  (itry  ,’eq  .  0)i  try  =  1 
m(i )=i try 
ib (i  )  =  volt0 
return 
end 

c _ 

c _ 

c 
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subroutine  sigavg 


c 

c 

c 

c 

c 

“outine  to  receive,  average  and  display  a  signal 
from  the  a-d  to  provide  a  trace  to  place  the 
cu  rsor s . 

integer*4  time1,time2 
integer  h r s , mi  ns , se cs , tck s 
byte  i chtab (2 ) 

real  array1(256),  array2(256) 

real  s x 1 ( 256 ) , sx2 ( 256 ) , sx2 1 ( 256 ) , sx22(256) 

common  /scalar/ic(2,4) , ib(2),m(2) 

common  /d at bu  f/i bu  fr ( 5 1 2 ), ibu  fr 1 (  256 ), i bu  fr2(  256 ) 

common  /stderr /se 1 ( 256 ) , se2 ( 256 ) 

common  /schmi t/i st , izzz 

common  ichtb 

equivalence (ichtb, ichtab) 

c 

c 

c 

c 

c 

variables:  sx1=  sum  of  xl 

sx21=  sum  of  x1**2 
sx2=  sum  of  x2 
sx22=  sum  of  x2**2 
call  p  lot 55 ( 9,0,0) 
call  plot55( 13,74, ) 
write (5, 1 000) 

1  000 

format ( ' $  of  signals  to  be  averaged:') 
accept  1500,  iternm 

1  500 
c 

format (i 1 0) 

2000 

write (5 , 2000) 

format ( ' $d o  you  wish  to  dump  the  buffer?  <y  or  n>:') 
accept  2500,idump 

2500 

c 

f  ormat (a  1  ) 

3000 

if  (ist.eq.1)  go  to  10  Iskip  if  schmidt  triggered 
wri te ( 5 , 3000) 

forma t (' $e nte r  the  time  between  samples  (sec):') 
accept  3500,  izzz 

3500 

c 

1  0 

format (i 5 ) 

do  20  i = 1 ,256 
array  1 ( i  )  =  0 
array2( i  )  =  0 
sx 1 ( i  )  =  0 
sx21 ( i  )  =  0 
sx2( i  )  =  0 
sx22 ( i  )  =  0 

20 

c 

c 

c 

continue 

collect  'iternm'  samples  and  a vg  them 

i count  =  256 
i ch  cnt  =2 
krt  =  3 
kcnt  =  1 

4000 

c 

c 

ich 1 =i chtab  ( 1 ) 
ich2=ichtab (2) 
do  90  itr=  1 , iternm 
call  gtim(timel) 

call  c o le ct ( i b ( 1 ) , i ch 1  )  !find  current  zero  level 

call  c ole ct ( i b ( 2 ) , i ch2 )  Ifind  current  zero  level 

if  (ist.eq.0)  call  stimul  ! i f  not  schmitt  trig,  run  trial 
call  getbu  f (i bu f r , i count , ichtab, ich cnt , krt , kcnt , ist , ierr ) 
write (5 , 4000)i tr 
format ( 35x ,  'work ing  ’,i4////) 

now  write  these  values  into  two  arrays  which  represent 
the  two  channels  of  input. 
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C 

C 

i  =1 
j  =  1 

ibu  fr ( 1 ) =i bufr (3  )  !first  point  is  garbage 

ibufr(2)=ibufr(4) 

30  if  (j.gt.256)  go  to  40 

reali =ibufr(i ) 

arrayl ( j )=array1 ( j )+(ibufr (i )-ib(1 ) ) 
sx  1  ( j  )  =sx 1 ( j )  +  ( reali -ib( 1 ) )/m( 1 ) 
sx21(j)=sx21(j)+(( reali -ib(1))/m(1))##2 

c 

reali =i bufr (i +1 ) 

array2( j )=array2( j )+(ibufr(i  +  1 ) -i b ( 2 ) ) 
sx2(j)=sx2(j)+(reali-ib(2))/m(2) 
sx22(j ) =sx22( j )+( (reali-ib(2) )/m(2) )  *  *2 
c 

J  =  j  +  1 
i  =  i  +2 
go  to  30 
c 

40  if  (ist.eq.1)  go  to  90  !skip  if  s.trig 

call  c  vtt imCtime 1 ,hrs,mins,secs,tcks) 

5-0  call  gtim(time2) 

call  cvttim(time2, ihrs , imins, isecs, itcks) 
if(isecs. eq . secs . and . itcks. eq.tcks)write(5,4500) 

4500  formatC*  turn  on  the  line  clock*) 

i f ( ( secs-i secs ) . It . i zzz . and . (i secs-secs ) . It . izzz )  go  to  50 
c 

c  go  get  next  sample 

90  continue 

ap  =i ternm 
do  60  i =1 ,256 

se 1 (i )  =  sqrt(abs(((sx21(i  )/ap)-((sx1(i )/ap)**2)/ap))) 
se2(i  )=sqrt(abs(((sx22(i  )/ap)-( (sx2(i )/ap)##2)/ap)))  !s. e. 

60  continue 

c 

i f (i dump . ne .  *y * )  go  to  80 
do  70  i  =  1  ,512,2 

70  write(5,5000)i,ibufr(i),ibufr(i+1) 

5000  f ormat ( i 5 , i 1 0 , i 1 0 ) 

write (5 ,5500 )m (1),ib(1),m(2),ib(2) 

5500  formate  m(  1  )r  * ,  i5, 3x  ,  *  i  b  ( 1 )  =  * ,  i5, // *  m(2)  =  * ,  i5,3x,  *ib(2)  =  »  ,15) 
stop 

c  plot  the  average 

80  call  axes 

call  trcplt (array  1 , array2  ,  iternm ) 
retu  rn 
end 
c 

c  _ 

c  . 

subroutine  trcplt (a r 1 , ar2 , itq ) 

c 

c 

c 

dimension  ar  1  ( 256 ) , ar 2( 256 ) 

common  /datbuf/ibufr (51 2 ) , ibufr 1 ( 256 ) , ibufr2( 256) 
common  /amplfy /i amp2 
common  ichtb,ilat 
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c 

ah=itq 

do  10  k=1  ,256 
ibufr 1 (k  )=ar 1 (k  )/ah 
ibufr2(k)=ar2(k)/ah 
1  0  continue 

20  call  plot 55 (9,45,0) 

call  plot55( 11,,) 

call  plot55(12,,’  amp  factor  :  ’  ) 

call  plot55( 13,72, ) 
do  15  i ai a= 1 , 59 
15  call  p lot55(  13,67,  ) 

accept  1000,  iampl 
1  000  format ( i 3  ) 

if  (iampl. eq. 0)return 
i amp2=i ampl 
c 

i 2=  1 0  .'signal  to  plot  on  lower  graph 

call  scaleCibufrl ,256,r1,r2,ilat) 
call  plot55( 1,0,) 
call  plot 55(7,1 11,10) 
call  plot (ibufr 1 ,256 , rl , r2, i2, iampl ) 
c 

i 2=  1 1 5  Isignal  to  plot  on  upper  graph 

call  scaled  bufr2, 256, rl  ,  r2 ,  i  la  t ) 

call  p lot55( 1,1,) 

call  plot 55(7 ,111,115) 

call  plot (ibufr 2, 256fr1,r2,i2, iampl ) 

go  to  20 

end 

c _ 

c 

c  routine  to  plot  the  data 

c 

subroutine  plot(i12,i1 , rl , r 2 , i2 , iampl ) 
c 

c  variable  names  consistent  with  'scale’ 

c 

integer  i 1  2 ( 25 6 ) 
real  r(256) 
if  ( r  1  . ne . r2 )  go  to  5 
write (5,1 000) 

1  000  format('  faulty  signal  read ’, /25x ,  ’terminate  execution’) 
stop 

5  r 1 = 1 00/ (r 1-r2)  !scale  data  to  fit  plot  size 

j  =  0 

do  10  i = 1 1 1 , 366 , iampl 
j  =  j+1 

r(  j )=i 1 2 ( j  ) 

10  call  plot55(8,i,  ifix(  (r  ( j)-r2)*rD+i2) 

return 
end 
c 
c 

c  ~  . 
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C 

c  this  routine  will  scale  the  data  arrays  for 

c  plotting,  it  returns  r1=max  value;  r2=min 

c  value,  r  is  the  input  array, 

c  il  is  the  of  points, 

c 

subroutine  scale (i 1 2 , i 1 ,r1 , r2f ilat ) 
integer  i 1 2 ( 256 ) 
dimension  r(256) 
do  10  i  =  1 ,256 
10  r  ( i  )  =  i  1 2  C  i  ) 

r 1  =  -1 . 0e  32 
r  2=  1 . 0e  32 
do  30  i =1 , i 1 
if  (r(i).gt.rl)  r1=r(i) 
if  (r(i).lt.r2)  go  to  20 
go  to  30 

2  0  r2  =  r ( i ) 

ilat=i  !  use  as  flag  for  location  of  stimulus  artifact 

30  continue 

retu  rn 
end 
c 
c 
c 

C  . . 

c  routine  to  place  cursors  on  the  traces  shown.  four 

c  cursors  can  be  placed  on  the  bottom  graph,  two 

c  cursors  are  placed  on  the  top  trace, 

c 

subroutine  csrplc 

common  /datbu  f /i bu  fr ( 5 1 2 ) , i bu  f r 1  ( 256 ) , i bu  fr2( 25  6 ) 
common  /stderr /se 1 ( 256 ) , se2 ( 256 ) 
c 

c  decide  which  graph  to  set  cursors  on. 

1 0  call  plot 55( 9,0,0) 

call  plot 55( 11,,) 

call  p lot 55 ( 12 , ,  * cu rsors  on  which  graph? (top  =  1 /bot  =  0)  2=exit:  ’ ) 
accept  1000,  igraph 

1  000  format (i 1 ) 

if  (igraph. eq .2)  go  to  70 
call  p lot 55 (9,0,0) 
call  plot 55( 1 1 ,, ) 
call  plot55( 1 , igraph,  ) 
c 

if  ( i graph. eq . 0 )  nc  =  4  ! i f  graphO,  get  4  cursors 

if  (igraph. eq. 1 )  nc  =  2  ! i f  graphl ,  get  2  cursors 

c 

c  identify  and  store  the  cursors  in  ’icf 

do  60  i = 1 , nc 

c  prompt  for  cursor  position 

20  call  p lot 55 (9,0,0) 

call  plot 55 ( 11,,) 

c  256=25.6  ms  ie:  crs  position=  10#  time 

call  plot55( 12 , , ’cursor  posi t i on ( 1 -256 )  use  <cr>  to  store:  f) 
call  plot55(9,43,0) 
accept  1500,  icp 
1 500  format (i 3 ) 

if  (icp. gt. 400)  go  to  40 
if  (icp. It. 0)  go  to  40 

if  (icp.eq.0)  call  st ore ( i , i graph , i last ) 
go  to  50 

40  call  p lot 55 ( 9, 0 , 0 ) 

call  plot 55(11 , ,  ) 

call  plot55( 12 ,, ’error  on  input,  hit  <cr>  to  continue’) 
accept  2000,  i i i 

2  000  format (a  1 ) 

go  to  20 

50  if  (icp.eq.0)  go  to  60 

call  plot55(6, ilast+1 1 0, 0) 
call  plot 55( 6 , icp  +  1 1 0 , 1  ) 

if  (i graph. eq . 0)call  d i sp ( i bu f r 1 , se 1 , i cp , i graph+1 ) 
if  ( i  grapji .  eq  .  1  )  ca  11  d i sp ( i bu fr2 , se2 , i cp , i graph* 1 ) 
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60 

70 

c 

c 

c 

G 

C 

C 

C 


c 

G 

C 

C 

C 

C 

C 

G 


1  000 


c 

c 

c 

c 

c 


c 


c 

c 

c 


ilast=iop 
go  to  2  0 
continue 
go  to  10 
return 
end 


this  routi he “enters  tHe  cursor  values  into  the 
table  'ic'  which  holds  the  values  for  later  use. 
as  well,  the  values  stored  are  displayed  on  the 
screen  as  well  as  the  cursor  positions. 

subroutine  store (i , i graph , ilast ) 
common  /scalar /ic(2, 4) , ib ( 2 ) , m(2 ) 
common  /amp lfy /i amp2 

i c ( ( i grap h  + 1 ) , i ) =i last /i amp  2  !the  last  value,  before  a 

zero  was  entered  to  exececute 
a  save,  is  placed  in  ic(i,j) 
where  *i  '  is  the  graph  and 
1  j  '  is  the  jth  cursor  stored. 

now  ,  display  the  saved  value  in  the  app.spot 

iplace=1-i graph 

call  plot 55(9,9, ((15*iplace)+2)+i ) 

wri te ( 5 , 1 000 )i last 

format (i 4 ) 

i last =0 

retu  rn 

end 


subroutine  axes 


call  p lot 55 ( 9,0,0) 
call  p lot 55 ( 10, , ) 

call  p lot5 5 C 1 2 f  f  *  note:  each  point  =  0.1  ms') 

call  plot55( 4, 1 ,115) 

call  plot55( 4 , 1 ,10) 

call  plot55(5, 370, 1 ) 

call  p lot 55 ( 5 , 1 1 0 , 1 ) 

call  plot55( 9 , 1 , 1 ) 

call  plot55( 12 ,, ’cursor  x') 

call  plot 55( 9,3,3) 

call  p lot 55 ( 12, ,  *  1 » ) 

call  plot 55( 9,3,4) 

call  plot55( 12, , ’2’ ) 

call  plot55( 9,1,16) 

call  plot  55 ( 1 2 , ,  ' cu rsor  x') 

call  plot55 ( 9,3,18) 

call  plot 55(  12, ,  ’ 1  '  ) 

call  plot 55(9,3,19) 

call  plot55( 12 , ,  ' 2'  ) 

call  plot55( 9,3,20) 

call  plot55( 12,  ,  *  3 ?  ) 

call  p lot 55 (9,3,21) 

call  plot55( 12 , ,  ’  4 '  ) 

call  plot55( 9 , 64 , 4  ) 

call  plot 55 ( 1 2 , ,  ' cu rrent  position') 

call  p lot 55( 9,64,5) 

call  plot 55( 12 , ,  'y =  *  ) 

call  plot 55( 9,64,6) 

call  plot55( 12, , 'se. y=' ) 

return 

end 
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C 

C 

subroutine  di sp (r , rx , iq r , iqx  ) 

dimension  r ( 256 ) , rx ( 256 ) 

common  /scalar/ic (2 , 4 ) , ib(2) , m(2) 

call  plot55( 9,68,5) 

write(5,1000)r(iqr)/m(iqx) 

call  plot 55(9,68,6) 

write ( 5 , lOOO)rx(iqr) 

1000  format (flO. 3) 
return 
end 
c 
c 


C . 3  SWAN RU .FOR  -  EXPERIMENT  RUN  SECTION 


c 

c 

c 

c 

c 

c 

c 

c 

c 


1  0 
20 


c 


950 

960 

1  000 
1  500 
25 


2000 

27 

2500 


3000 

3500 

4000 


subroutine  call  mapping: 

main  program 

/ 

swanset  runit 

/ 

graph  runex 

/  i 

getbuf  tvmax  negmax 

subroutine  runit 
byte  i f i le ( 1 4 ) 
byte  blank 

common  /scalar/i  c (2, 4 ) , ib ( 2) , m(2) 
common  /u ni ts/iuni ts , iu ni t2 

common  /outf il/laten (  1 000 ) , mvolt (  1000) , negex 1  (  1000) ,  nege x2(  1 000 ) 

common  ichtb, ilat , idi f , itest 

data  blank/’  '/ 

do  20  i  =  1 ,  1  4 

ifi le (i )=blank 

call  plot55(2, 1+2+4+32+64+51 2,) 

call  plot 55 (9,0,0) 
call  p lot 55 ( 10,,) 
write(5,950) 

formatC 1 $d i ff erences  will  be  relative  to  zero(0)  or  1st  cursor(l):  ') 
accept  960,  idif, itest 
format (2i2 ) 
write (5,1 000 ) 

formatC1  supply  the  information  requested  on  the  max.  values') 
format(///,20x,  ’ post synap t i c ' , //,  ’$maxneg  cursors  1-2  (  ,  [mv/uv])  :  f) 
write(5, 1500) 
accept  2000 , isl , iuni ts 

i f (iuni ts . ne .  'uv ' . and . iuni ts . ne.  Tmv * ) go  to  25 

if (iunits.eq.  'uv' )iunits  =  1000  fscale  change  to  microvlt 

if(iunits.eq.  'mv')iunits  =  1 

format (i 6 , a2 ) 

write (5 ,2500) 

format (//, 20x ,  'presynaptic ',//,  1  maxneg  cursors  1-2  (  .[mv/uv]  :  '.$) 

accept  2000 , is2 , iuni t2 

if ( iuni t2  .  ne .  'u v '. and . iuni t2  .  ne .  'mv  '  )go  to  27 
i f ( iuni t2 . eq . 'u v * ) i uni t2= 1 000 
if(iunit2.eq.  *mv')iunit2=1 
f  ormat ( i 5 ) 
wri te ( 5 , 3500 ) 

format(//,'  maxneg  ex (mv  ); cursors  3-4:  ',$) 

accept  2000, is3 
write ( 5 , 4000 ) 

format (/,' $how  many  trials?  (  1000  max):') 

accept  3000,iiiii 

call  graph (i s 1 , i s2 , i s3  ) 

call  runex  (isl , is2 , is3  ,  iiiii  ) 
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c  write  the  results  for  today  into  a  storage  file 

call  p  lot55( 9,0,0) 
write (5,^500) 

4500  format (5x ,$ 'enter  file  name  for  storage-  dkl  :<date>.dat  :') 
accept  5000,ifile 

5  000  f ormat ( 1 4a 1 ) 

open  (uni  t  =  1 ,  name  =  i  f  ile,type='new') 
do  30  i  =  1  ,  iiiii 

aa  =  f  loat  (mvolt  (i  )  ) /float  (m_(2)  ) 
ab=float(negex1(i ) ) /float (m ( 1 ) ) 
ac=float (negex2(i ) ) /float  Cm  C 1 )) 
write(1,5500)laten(i ), aa , ab , ac 
30  continue 

5500  format (i 5 , 3f  1  0. 4) 

i f lag=-5 

write (1 ,5500)iflag, aa,  ab, ac 
close(unit  =  1 ,dispose='save  * ) 
call  p lot 55( 9,0,0) 
call  plot 55( 10, , ) 
write ( 5 ,6000) 

6000  format ( ’ $do  you  wish  to  run  more  trials  with  the  same  setup?') 
accept  6500, iask 

6  500  f ormat (a  1 ) 

if (iask . eq . 'y 1 )  go  to  10 

return 

end 


c  this  routine  sets  up  the  graphics  display  for 

c  the  experiment 

c 

subroutine  graph (i s 1 , i s2 , is3 ) 

common  /scalar/i c(2, 4 ) , ib(2) , m(2) 

common  /uni ts/iuni ts , iuni t2 

common  i chtb, ilat 

call  plot 55( 9,0,0) 

call  plot55( 10, , ) 

call  p lot 55 ( 9,20,0) 

call  plot 5 5 ( 12 ,,  'total  elapsed  time:') 

call  plot 55( 4  , 1 ,30) 

call  plot 55 (  4,1 ,80) 

call  plot 55( 4 , 1 ,129) 

call  plot 55( 4 , 1 ,1  30) 

call  plot55( 4 , 1 ,180) 

call  plot 55 (5,111,1) 

do  20  j  =  5 , 20 , 5 

do  10  i = 1 8 , 79 , 3 

call  plot55( 9 , i , j) 

1  0  call  plot55( 12, ,  ’  +  ' ) 

20  continue 

do  30  i =5 , 20 , 5 
call  plot55( 9,74, i ) 
call  plot 55( 12, , 'trials ’ ) 
call  p lot 55( 9 , 47 , i  ) 

30  call  p lot 55 ( 12, ,  ' 100'  ) 

call  plot55(9, 12,2) 

write(5,1000)(ic(2,2)-ilat)/9  !1/2  of  the  scaled  values 

call  plot 55( 9,12,6) 
write ( 5 , 1000)is1/2 
1 000  format(i6, ' ' ) 

call  plot55( 9, 12,11) 

write ( 5 , 1000)is2/2 

call  plot 55 (9 ; 12 , 16 ) 

write (5,1 000)is3/2 

call  plot 55 (9,0,0) 

call  plot 55( 1 2 , , 'p  latency') 

call  plot55( 9,0,1  ) 

call  plot55(  12,  ,  'o  stimulus->max. ' ) 

call  plot 55(9,0,2) 

call  plot 55 ( 1 2 , , ' s  neg.excurs') 
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C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


call  plot55(9,0,3) 

call  p lot55( 12, , ’t  (msec) ') 

call  plot 55 (9,0,4) 

call  plot55( 12 , , ’s ’ ) 

call  plot 55 (9,0,5) 

call  plot 55 ( 1 2 , ,  ’y  delta  voltage’) 

call  p lot 55 (9,0,6) 

call  plot  55(  12 ,  ,.’n  max.’) 

call  plot55(9,0,7) 

call  plot55( 12, , ’p  neg.excurs’) 

call  plot 55(9,0,8) 

if (iuni ts . eq . 1 )call  plot 55 ( 12 , ,  ’s  (mv)') 

if (iunits. ne. 1 )call  plot 5 5 ( 12 , ,  's  (uv)') 

call  plot 5 5 ( 9,0,11) 

call  plot 55 ( 12 , ,  ’p  delta  voltage’) 

call  plot 5 5 ( 9 ,0,12) 

call  plot 55( 12 , , *r  crsrs  1->2’) 

call  plot55( 9,0,13) 

i f (iuni t2 . eq . 1 )eall  plot55 (  12  ,  , ’e  (mv)') 

i f ( iuni t2 . ne. 1 )ca 11  plot 55 ( 1 2 , , ’e  (uv)’) 

call  p lot 55 ( 9,0,16) 

call  plot55(  12, ,  ’s  delta  voltage') 

call  plot 55 (9,0,17) 

call  plot 55( 12 , ,  *y  crsrs  3->4’) 

call  p lot 55 ( 9,0,18) 

call  plot55( 12 , ,  fn  (mv)') 

call  plot 55 (9,0,19) 

call  plot55( 12, ,  ’p  ’  ) 

call  plot 55( 9 , 0 ,20) 

call  plot 5 5 (  12 , ,  ’s’) 

call  plot55( 9,42,22) 

call  plot 55( 12, ,  ’trial:  '  ) 

return 

end 


routine  to  run  the  experiments  and  to  collect 
and  display  the  data  on  the  graphs 

subroutine  runex(is1 ,is2,is3,iiiii) 

byte  ichtab(2) 

common  /d  atbu  f /i  bu  fr ( 5 1 2 ) , i bu  fr 1  ( 256 ) , ibu  fr2( 256 ) 

common  /scalar/ic(2,4), ib(2) , m(2) 

common  /howmny/i ter ,  item 

common  /schmit/ist , izzz 

common  ichtb,ilat 

integer*4  time1,time2 

equi valence (ichtb, ichtab) 


variables: 

i  c 


the  table  of  cursor  values 
ib,m  the  parameters  for  output  transform 
iter  trial  counter 

item  trial  counter  for  display  control 

laten  array  for  storing  latency  values  , ch  2 

mvolt  array  for  storing  max.  neg.  values, ch  2 

negexl  array  for  storing  max.  neg.  values, 1st  curs. 

negex2  array  for  storing  max.  neg.  values, 2nd  curs, 

ibufrl  data  buffer,  ch  1 
ibufr2  data  buffer,  ch  2 

time1,time2  current  time  for  sampling  control 

izzz  time  between  samples  (snoozzzzz  time) 


,  ch 
,  ch 


1 

1 
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C 

iter=0  ! set  initial  value 
item  =  1 
icou  nt  =256 
i chcnt =2 
krate=3 
kcount  =  1 
ichl =ichtab  ( 1 ) 
ich2  =  i chtab (2 ) 
c 

10  call  gt im ( time  1 ) 

call  cvttimCtimel ,ihrs,imin,isec1,itck) 
c 

c  run  clock  at  lOkhz,  convert  at  ea .  overflow 

if(ist.eq.O)  call  stimul  !run  the  trial  if  not  schmitt  trig, 
call  colect (ib( 1 ), ichl )  !get  current  zero  level 

call  colect (ib (2) , ich2 )  ! get  current  zero  level 

call  get bu f (ibufr, icount , ichtab,ichcnt,krate, kcount, ist, ierr ) 
i  =1 

J  =  1 

ibu  f  r  ( 1  )  =i  bu  fr  (  3 ) 

i bu f r ( 2 ) =i bu f r ( 4 )  !  first  points  are  junk 

18  if (i . gt . 256)go  to  20 

ibuf rl  (i  )=ibufr  ( j  ) 

ibu f r2( i ) =i bu f r ( j+1 )  !separate  both  channels 

j=j+2 
i  =i  +1 
go  to  18 

20  iter=iter+1 

call  p lot55( 9 ,50,21) 
write(5, 1 000)iter 
1  000  f  ormat ( i 5 ) 

c 

if  (iter. eq . iiiii  )  go  to  50 
if  (iter. eq . 199)  call  copy 
if  (i ter. eq . 399)  call  copy 
if  (i ter . eq . 599 )  call  copy 
if  (i ter. eq . 799 )  call  copy 
call  plot55( 5, itern+1 1 1 , 0) 
call  plot55(5, item +  113,1) 
c 

call  tvmax (ibu fr2 , isl )  Iprocess  ch  2 
call  negmax (i bufr 1 , is2 , is3 )  !process  ch  1 

call  p lot 55( 9,0,23) 
call  p lot 55( 10, , ) 
item  =i tern+2 
c 

c  wait  for  izzz  seconds  to  elapse  before  running  the  next  trial 

30  call  gtim(time2) 

call  cuttim(time2,ihrs, imin,isec2,itck) 

if(ist.eq.l)  go  to  40  !skip  the  wait  if  schmitt  triggered 
if  ( (i sec2-i sec  1 ) . It . izzz. and . (isec 1-i sec2) . It.  izzz )  go  to  30 
40  call  plot55( 9,41 ,0) 

writae(5,1500)ihrs, imin, isec2 
1  500  forma t (i2 ,  ,:t, i2 i2  ) 
go  to  10 
c 

50  call  copy 

return 
end 

- - 

c 

subroutine  copy 
common  /howmny /i ter,  item 
call  plot55(9,79,23) 
call  plot 55( 13,93, ) 
i tern  =  1 
return 
end 
c 
c 

c  ZZZZ  _ ; — - - .  - - 
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C 

c  this  routine  uses  the  cursors  designated  for  the  post- 

c  synaptic  side  to  find  the  most  negative  value  sampled 

c  between  the  cursors  and  the  time  when  this  value  occured. 

c  this  time  is  measured  from  the  initial  depol.  and  is 

c  displayed  as  a  latency, 

c 

subroutine  tvmax (i store , isl ) 

dimension  istore(256) 

common  /scalar/ic(2,4),ib(2),m(2) 

common  /outfi l/laten(1000), mvolt (1000), negex1(1000),negex2(1000) 
common  /howmny /i ter,  item 
common  /units/iuni ts , iuni t2 
common  i cht b , ilat , id i f , i t est 
i vmi n=1 0000 
c 

do  10  i=ic(2, 1 ) , ic(2,2)  ! look  between  cursors 
if(itest.eq. 1 ) wr i te ( 5 , 1000)istore(i  ) 

1  000  format (i 1 0 ) 

if  (istoreCi ). gt.ivmin)  go  to  10 

i vmi n  =  i store ( i  )  difference  relative  to  zero 
lat  =i 

1  0  conti nue 

if (itest . eq . 1 )wri te ( 5 , 1 500 )i vmi n , lat 
1500  formatC*  minimum= 1 , i 1 0 ,  *  at  lat  =’,i10) 

i f (idi f . eq . 1 )i vmi n=i store (i c (2 , 1 ) )-i store (lat )  !diff  rel  to  1st  curs 
if (i vmin . It . 0)i vmi n=-i vmin 
c 

laten (i ter )=f loat (lat-i lat )/ 4 .65 

if ( (lat . ne . ic (2, 1 ) ). and. (lat . ne. ic (2,2) ))go  to  20 
la t =i  la t 
i vmi n  =  0 

20  mvolt (iter  )  =  i vmin 

c  graph  it  now 

i y 1 =( (50 . / (f loat (ic(2,2 )-i lat )/4.65) ) * ( float ( lat-i lat ) / 4 . 65 ) ) + 1 80 . 

if  (iyl . gt . 230 )iy 1 =230 

if  (iyl . It . 1 8 0 ) i y 1  =  1 80  !  check  scaling 

c 

iy2= ( 50 . /float (isl) ) *i u ni ts * ( float (i vmi n)/float(m(2)))+130. 

if  (iy2 . gt . 1 80)iy2= 1 80 

if  (iy2.lt . 130)iy2=130 

call  p lot 55 ( 1 , 1 ,  ) 

call  plot55( 3, 1 1 0+itern , iyl  ) 

call  plot55( 3 i 1 1 1+i tern , iy2) 
return 
end 
c 

c  - - - - - - - - - - — 

c  - - - — - - - - - - - —  - - 

c 

c  this  routine  finds  the  most  neg.  points  within  the  time 

c  windows  provided  by  cursor  sets  1-2  3-4.  the  input 

c  examined  is  from  channel  1,  presynaptic. 

c 

subroutine  negmax (i store , is2 , is3) 

dimension  istore(256) 

common  /scalar /i c (2 , 4 ) , ib (2 ) , m(2) 

common  /outfil/laten( 1000), mvolt ( 1000) , negexl ( 1000) , negex2( 1000) 
common  /howmny/i t er ,  item 
common  /uni ts/iunits , iuni t2 
common  i chtb, ilat , idif , itest 
c 

do  20  j  =  1  ,  3 ,2 
i vmi n  =  10000 

do  10  i =i c ( 1 , j ) , ic ( 1 , j+1 ) 
if(itest.eq. 1 )wri te( 5 , 1 000 )i store (i ) 

1000  format(ilO) 

if  (i store (i  ). ge. i vmi n )  go  to  10 
i vmi n  =i store ( i  ) 
lat  =i 
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1 0  continue 

if(itest.eq.1)write(5, 1500 )i vmi n,  lat 
1  500  format(?  minimum  =  f , i 1 0 , f  lat=’,i10) 

ifCidif.eq.l )i vmin=i store (ic ( 1 , j ) )-i store (lat ) 
i f (i vmi n. It . 0)i vmi n=-i vmi n 
c 

if  (j.lt.2)  ne gex 1 ( i ter ) =i vmi n 
if  (j.gt.2)  negex2( i ter  )  =i vmi n 
call  plot 55( 1,0,) 
if  (j.gt.2)  go  to  30 
c 

if (lat. eq . ic( 1 , j+1 ) . or . lat . eq . ic (1 , j) )  i  vmin  =  0 
c  !if  the  min  pt.@  curs.,  defined  as  0 

iz 1=( 50. /float (is2)*iunit2)#( (float (i vmi n ) /float (m ( 1 ) ) ) )+80. 
if  (i z 1 .gt . 130)iz1  =  130 
if (izl . It . 80 )i zl  =  80 
call  plot55( 3, 1 10+i tern , izl ) 
go  to  20 
c 

30  iz2= ( 50 . /float (is3) )*( (float (i vmi n ) /float (m ( 1 ) ) ) ) +30 . 

if (iz2.gt.80)iz2=80 
if (iz2.lt. 30)iz2=30 
call  plot55( 3, 1 1 1+i tern , iz2) 
c 

20  continue 

return 
end 
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D.1  COLECT  .  MAC 


title  colect 
ident  /donOI/ 


d.o.  norman 


routine  to  take  a  single  a/d  on  a 
given  channel 


fortran  calling  procedure 

call  colect (dataword , channel ) 

dataword  :::  receives  the  result  of  the  a/d 

channel  : : :  input  channel 


colect 
s  tart : 


wait: 


adsr  = 

170400 

adbr  = 

170402 

clr 

@  adsr 

; clear  a/d  stat  reg 

mo  v 

2(r5) , rO 

;set  up  address  for  storage 

mo  v 

@4(r5),r1 

;get  channel  number 

a  si 

r  1 

asl 

r  1 

shift  it 

asl 

r  1 

into  the 

asl 

r  1 

correct 

asl 

r  1 

position 

asl 

r  1 

before 

asl 

r  1 

’ ori ng  9  it 

asl 

r  1 

with  the  adsr 

bi  s 

r  1 , @  adsr 

set  up  channel 

inc 

@  adsr  ; 

;  start  conversion 

ts  t  b 

@  adsr  ; 

;  done? 

bpl 

wait  ; 

if  not  then  wait 

mo  v 

@  adbr, (rO)  j 

;  load  data  word 

rt  s 

pc 

.  end 

D  .  2  ST IMUL  .  MAC 


title  stimul 
ident  /don02/ 
d.o.  norman 

routine  to  cause  a  triggering  of 
the  prep  stimulator  through  the 
external  triggering  connection. 

assumes  a  drv-11  (parallel  output) 
is  used  and  set  up  according  to 
specs . 

trigger  is  through  bit  0 

fortran  call: 
call  stimul 
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stimuli : 


dr vou  t  = 

167772 

clr 

@  drvout 

inc 

@  drvout 

nop 

nop 

nop 

nop 

nop 

nop 

nop 

nop 

nop 

nop 

clr 

@  drvout 

rt  s 

pc 

.  end 

;  toggle  bit  0  (turn  it  on) 
;  wait  a  while 


;  clear  bit  0  (turn  it  off) 


.title  getbuf  -  rtll  adv/kwv  driver  module  -  program  control 
; ;  .sbttl  getbfl  -  burst  mode  sampling  version 

. ident  /lmf 01  a/ 

.nlist  end 
.enabl  lc 

. mcall 

written: 

1.  m.  fraser  21-mar-80 
call  syntax:  ...fortran 


call 

get  bu  f  (i  bu  f ,  icount ,  ichtab,  ichcnt,krate,kcount,  iwait ,  ierr) 


where : 

ibuf 

icou  nt 
ient ab 
interest . 

ichcn t 
kra  t  e 

manual 

kcou  nt 

iwait 

ierr 


integer  buffer  to  receive  data  in  order  as  sampled. 

must  be  dimensioned  i count *i chcnt  long 
number  of  sample  sets  to  convert 
*byte*  table  of  channel  numbers  in  order  of 

length  of  channel  list 

rtc  clock  count  units  parameter  (1  =  mhz)  see  users 

number  of  units  above  between  overflows 
wait  for  st2  trigger  if  ><  0 

error  indicator  (reflects  only  data  overruns) 


define  parameter  below  to  select  level  of  error  indication 

errmrk  not  defined:  no  error  checking 

errmrk  =  0  count  overflows  only 

errmrk  =  1  count  errors  set  data  point  bit  12  on  error 

errmrk  =  1 
local  macros 

parameter  error  reporter  -  default  error  81.  "invalid 
argument" 


.macro  prmerr  prm 
mo  v  prm, prm. 8 ( r 5 ) 

mo  v  spsa  ve , sp 

return  ;  return  with  error  report 

.endm  prmerr 
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get  address  from  fortran  argument  list  or  check  of  loc 
blank  (ok  if  ne) 


.macro 

geta  off, loc, ?a , ?b 

.iif  idn  loc,r5  .error  ; you ’ 

(r5) 

cmpb 

off/2, (r5) 

ble 

a 

prmerr 

off 

a  : 

V 

.if  nb 

loc 

mo  v 

off  (r5) , loc 

cmp 

loc,  -1 

bne 

b 

trap 

128. +16. 

.  iff 

cmp 

of f (r5 ) ,  -1 

.  endc 

b  : 

.  endm 

geta 

> 

;  get 

parameter 

from  fortran  argument 

only 

> 

. macro 

getp  off, loc 

geta 

of  f , loc 

mo  v 

§loc , loc 

.  endm 

getp 

;  put 

parameter 

to  fortran  parameter 

. macro 

pu tp  off, loc , ?a , ?b 

cmpb 

off/2, (r5) 

ble 

b 

prmerr 

off 

b: 

cmp 

off (r5) ,  -1 

bne 

a 

9 

prmerr 

off 

Q  « 

.if  nb 

1  oc 

mo  v 

loc, off (r5) 

.iff 

clr 

off (r5) 

.  endc 

.  endm 

.  page 

;  local  definitions . 

word  arguments 


define  fortran  parameter  offsets 


prtn.1 

= 

2 

p  rm.  2 

= 

4 

p rm.  3 

= 

6 

p rm.  4 

= 

10 

p  rm.  5 

= 

12 

p rm.  6 

= 

14 

prm.  7 

= 

16 

p  rm.  8 

= 

20 

;  define 

ad v 1 1  and  kwv 1 1 

a  dcsr 

r 

170400 

a  sbu  f 

= 

adcsr  +  2 

k  wcsr 

= 

170420 

k  wbpr 

= 

kwcsr  +  2 
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p  r  0  =0 

pr7  =  340 


define  adv  csr 


ad .  go 

- 

1 

a  d . mad 

= 

4 

a  d . ide 

= 

10 

a  d . exe 

= 

20 

a  d . cle 

= 

40 

a  d . ine 

- 

100 

a  d . don 

= 

200 

a  d . adr 

= 

400 

a  d . ere 

= 

40000 

a  d . err 

= 

100000 

i 

a  d . idn 

= 

10000  ; 

;  define  kwv  csr 

k  w .  go 

— 

1 

kw . mod 

= 

2 

kw . rat 

= 

10 

k w . ine 

= 

100 

k w . ovf 

= 

200 

k  m.  st  1 

= 

400 

k  m . st2 

= 

1000 

km . osc 

= 

2000 

k  m. di o 

= 

4000 

k w . ovr 

= 

10000 

k w . gs2 

- 

20000 

kw. s2e 

= 

40000 

k  w . st2 
*  page 

= 

100000 

;  local 

s  torage 

bit  offsets 
start  immediately 

maint  -  load  all  data  fr#om  mux  bit  0 
set  identify  mode 
external  start  enable 
clock  overflow  start  enable 
enable  interupt  on  done 
done/ready/i nt erupt  request 
base  of  mux  address  (high  byte) 
enable  interrupt  on  error 
error  interrupt  request 

data  buffer  identify  bit 

bit  offsets 

start  counting  immediately 
base  bit  of  count  mode 
base  bit  of  rate  specifier 
interupt  on  overflow  enable 
overflag  flag/done 
maint  -  simulate  stl 

maint  -  simulate  st2 

maint  -  simulate  oscillator  single  cycle 

maint  -  disable  internal  oscillator 

error  overrun  flag 

enable  start  on  st2  overflow 

interrupt  on  st2  enable 

st2  set  flag/  interrupt  request 


s  psa  ve : 

.  word 

0  ; 

saved  sp  for  error  return 

c  ht ab : 

.  word 

0  ; 

storage 

for  channel  table  address 

c h cnt : 

.  word 

0 

temp  stoage  for  channel  count 

k  wfla  g: 

.  word 

0  ; 

rate  mode  words  for  clock 

dump : 

.  word 

0  ; 

place  to  dump  junk 

;  entry 

poi  nt 

g  e  t  bu  f : 

clr 

§  adcsr 

clr 

§  kwcsr 

mov 

sp  ,  sp  sa  ve 

save  sp  for  error  return  on  this 

geta 

prm. 3 , r2 

get  address  of  channel  table 

mov 

r2,  rO 

copy  addr  for  later 

mov 

r2, chtab 

save  address  again 

getp 

prm. 4 , r3 

get  number  of  channels 

mo  v 

r  3,  r  1 

save  for  later 

mov 

r  3 , chcnt 

save  count  again 

check 

all  channels  for 

validity 

1  $: 

bi  t  b 

C17, (rO) 

+ 

check  for  any  high  bits 

beq 

10$ 

continue 

sub 

r 2 ,  rO 

calc  element  in  error 

add 

10000, rO 

indicate  not  a  statement  no 

prmerr 

rO 

force  error  call  traceback 

1  0$: 

dec 

r  1 

check  all  of  list 

bne 

U 
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geta 

prm. 1 , rO 

get  address  of  data  buffer 

getp 

prm. 2 , r 1 

get  count  of  samples 

mov 

r 0 ,  r4 

copy  address 

ror 

r  *4 

get  low  bit 

bcc 

15$ 

odd  address  if  set 

prmerr 

prm.  1 

force  bad  buffer  address  error 

5$: 

getp 

prm.  5  ,  r*4  ; 

;  check  clock  rate 

bi  t 

C7,r4  ; 

;  check  clock  rate  (only  1-7) 

beq 

20$  ; 

;  ok  if  equal 

prmerr 

prm. 5  ; 

;  flag  parameter  number 

1 0$ : 

asl 

r  *4 

shift 

asl 

r  *4 

into 

asl 

r  *4 

position 

mov 

r*4 ,  kwf  lag 

save  into  temp 

getp 

prm.  6  ,  r *4 

get  number  of  intervals 

neg 

r  *4 

clock  counts  up  to  zero 

mo  v 

r4,@  kwbpr 

load  clock  preset  buffer 

bis 

kw. mod  *  1 , kwf lag 

;  set  to  repeat  interval 

mo  v 

kw. go ,  —  ( s p  ) 

;  assume  no  wait 

getp 

prm.7,r*4 

;  get  flags  word 

bne 

30$ 

;  good  guess 

mo  v 

kw. gs2 , (sp  ) 

;  set  up  for  wait  for  st2 

;  0$ : 

bis 

(sp )  +  , kwflag 

;  set  into  mask  word 

assume  no  conversion  errors 

putp 

prm. 8  ; 

;  clear  error  indicator 

mov 

ad . cle , §  adcsr  ; 

;  enable  clock  starts 

page 


note : 


;  the  following  code  runs  at  priority  7  (  non-interrupt ible ) 

;  to  make  data  acquisition  as  fast  as  possible.  to  drive 
the  advll-aa 

;  at  its  fastest  rate  the  loop  below  must  complete  within 
abou  t 

;  *40  microsec  nominal.  this  includes  9  ms  settling  time  for 

the  amps  and 

;  32  ms  conversion  time.  at  this  point  all  parameters  have 

been  validated 

;  and  set  up  in  the  following  registers: 

;  rO  -  current  buffer  pointer  (next  data  goes  here) 

;  rl  -  number  of  samples  remaining 

;  r2  -  pointer  to  channel  table  (next  channel  to  be  sampled) 

;  r3  -  channels  done  in  this  scan 

;  r*4  -  adc  csr  address 

;  r5  -  fortran  argument  pointer 


main  sampling  loop... 


mtps  pr7  ;  lock  out  interrupts 

mov  adcsr,  r *4 

mov  kwflag,@  kwcsr  ;  ok,  start  up  clock 


b loop : 

mo  vb 

(r2)+, 1 (r4) 

1  oop : 

tstb 

§  adcsr 

bpl 

loop 

mov 

2(r *4) ,  dump 

clr 

2(  r  *4 ) 

i  nc 

§  adcsr 

;  ;  set  up  channel 
;;  wait  for  adc  to  finish 
; ;  loop  until  done 
; ;  dump  fi rst  con v. 


; ;  then  get  another 
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1 ooop : 


1  0$: 


20$: 


tstb 
bpl 
mo  v 

@  adcsr 
looop 

2(r4), (r0)+ 

J  S 

;  ;  save  away  data 

mo  v 

2( r4) , dump 

> 

:  ;  ensure  bit  7  unset 

clr 

dec 

2  ( r  4 ) 
r3 

t 

;  count  down  channels 

beq 

10$ 

j 

;  continue 

mo  vb 

(r2) +, 1 ( r4) 

;  set  up  next  channel 

i  nc 

@  adcsr 

> 

;  start  up  adc 

br 

loop 

;  go  again 

rrmrk 

;  ;  ;  errmrk  - 

flag 

data  overruns  if  eq  0, 

mark 

eq  1 

.if  eq 
bi  c 

errmrk- 1 

ad . ide, @r 4 

>  i 

;  ;  assume  no  error 

.  endc 
tst 

@r  4 

J  ! 

;  ;  data  overrun?? 

bpl 

20$ 

»  ! 

;  ;  nope 

.if  eq 
bi  s 

errmrk- 1 
ad . ide , @r 4 

J 

;;  set  to  mark  data  in 

error 

.  endc 
dec 

@prm. 8(r5) 

) 

;  ;  count  errors 

.  endc 

; ; ;  errmrk  - 

flag 

e  rrors 

mo  v 

chtab, r2 

> 

;  reset  channel  table 

addre. 

mo  v 

chcnt , r3 

;  reset  channel  count 

dec 

H 

> 

;  count  down  sample  sets 

bne 

bloop 

> 

;  go  again 

)ne  now 

clr 

§  adcsr 

; ;  reset  adc 

clr 

@  kwcsr 

;  ;  and  clock 

mtps 

prO 

allow  interrupts 

retu  rn 
.  end 


;  return  to  main  program 
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APPENDIX  E.  ARCHITECTURE  OF  DATA  STORAGE  FILE 


To  allow  the  data  to  be  viewed  by  the  investigator  using  the  available  utility  routines 
(pip,  etc.),  the  data  are  stored  as  a  sequence  of  four  entries  per  trial.  As  currently 
implemented,  the  first  entry  is  the  latency  to  the  most  negative  excursion  within  the 
cursors  in  the  first  channel,  and  the  second  is  the  magnitude  of  the  excursion  in  mV. 
The  third  and  fourth  entries  refer  to  the  most  negative  excursions  between  the  cursor 
sets  on  the  second  channel.  The  "end  of  file"  mark  is  implemented  as  a  negative 
latency. 

The  data  are  written  into  a  file  with  "formatted"  and  "sequential  access"  attributes 
as  a  series  of  I5,3F10.4  records. 
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APPENDIX  F.  PROGRAMS  FOR  AUXILIARY  ANALYSIS 


At  present,  one  auxiliary  program  exists  that  uses  the  data  files  created  to  perform 
subsequent  analyses  of  the  data. 

PROGRAM:  REG.FOR 

PURPOSE:  To  determine  interactions  between  variables 


METHOD:  Scatterplot  of  variables  shown  on  VT55  screen  and  display  of  Pearson's  r 

calculated  for  the  data  set 
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